Intra-day rolling scheduling method for integrated heat and electricity system

ABSTRACT

An intra-day rolling scheduling method for an integrated heat and electricity system is provided. The method includes: establishing an objective function for scheduling of the integrated heat and electricity system, the objective function aiming to make operating costs of the integrated heat and electricity system to be a minimum; establishing constraints for a steady-state safe operation of the integrated heat and electricity system; and solving the objective function based on the constraints by an interior point method, to obtain an active power and a heating power of each combined heat and power unit, an active power of each thermal power unit, a heating power of each heat pump, and an active power consumed by each circulating pump, as an intra-day rolling scheduling scheme of the integrated heat and electricity system.

CROSS-REFERENCE TO RELATED APPLICATION

This is a continuation of International Application No. PCT/CN2017/114465 filed Dec. 4, 2017, which claims priority to Chinese Patent Application No. 201710989015.9, filed Oct. 22, 2017, the entire disclosures of which are incorporated herein by reference.

FIELD

The present disclosure relates to an intra-day rolling scheduling method for an integrated heat and electricity system, belongs to the field of operation and control technologies for a power grid with a plurality of energy forms.

BACKGROUND

The multi-energy system may include a plurality of energy flow subsystems. The energy may refer to such as electricity, heat, cold, gas, or transportation. The multi-energy system has considerable benefits, such as reducing energy consumption, improving overall energy efficiency, and reducing total energy costs, but improves the control complex.

The integrated heat and electricity system is one of the multi-energy systems. The integrated heat and electricity system may include a power grid and a heating network, in which coupling components such as combined heat and power units, heat pumps, and circulating pumps may couple the power grid and the heating network. However, there coupling components bring about new requirements on operation and control on the integrated heat and electricity system, such as the intra-day rolling scheduling on the integrated heat and electricity system needs to consider these coupling components comprehensively.

SUMMARY

The intra-day rolling scheduling method for the integrated heat and electricity system, provided by the present disclosure, may include the following: establishing an objective function for scheduling of the integrated heat and electricity system, the objective function aiming to make operating costs of the integrated heat and electricity system to be a minimum, the integrated heat and electricity system comprising combined heat and power units, thermal power units, heat pumps, and circulating pump; establishing constraints for a steady-state safe operation of the integrated heat and electricity system; and solving the objective function based on the constraints by an interior point method, to obtain an active power and a heating power of each combined heat and power unit, an active power of each thermal power unit, a heating power of each heat pump, and an active power consumed by each circulating pump, as an intra-day rolling scheduling scheme of the integrated heat and electricity system.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram of an integrated heat and electricity system according to an embodiment of the present disclosure.

FIG. 2 is a flowchart of an intra-day rolling scheduling method for an integrated heat and electricity system according to an embodiment of the present disclosure.

DETAILED DESCRIPTION

The integrated heat and electricity system is shown in FIG. 1. The integrated heat and electricity system may include a power grid and a heating network. Coupling components such as combined heat and power units, heat pumps, and circulating pumps may couple the power grid and the heating network.

The intra-day rolling scheduling method for the integrated heat and electricity system, provided by the present disclosure, may include: establishing (110) an objective function for scheduling of the integrated heat and electricity system, the objective function aiming to make operating costs of the integrated heat and electricity system to be a minimum, the integrated heat and electricity system comprising combined heat and power units, thermal power units, heat pumps, and circulating pump; establishing (120) constraints for a steady-state safe operation of the integrated heat and electricity system; and solving (130) the objective function based on the constraints by an interior point method, to obtain an active power and a heating power of each combined heat and power unit, an active power of each thermal power unit, a heating power of each heat pump, and an active power consumed by each circulating pump, as an intra-day rolling scheduling scheme of the integrated heat and electricity system.

In detail, the intra-day rolling scheduling method for the integrated heat and electricity system, provided by the present disclosure, may include the following steps.

(1) An objective function for optimal scheduling of the integrated heat and electricity system may be established by a formula of:

${\min {\sum\limits_{t = 1}^{16}{\left( {{\sum\limits_{b = 1}^{N}{F\ \left( {p_{b,t},\ q_{b,t}} \right)}} + {\sum\limits_{x = 1}^{N_{TU}}{F_{TU}\left( p_{x,t} \right)}}} \right)\Delta t}}},$

where, p_(b,t) represents an active power of a b^(th) combined heat and power unit in the integrated heat and electricity system at a t^(th) scheduling period, q_(b,t) represents a heating power of the b^(th) combined heat and power unit in the integrated heat and electricity system at the t^(th) scheduling period, N represents a number of combined heat and power units in the integrated heat and electricity system, F(P_(b,t),q_(b,t)) represents a running cost of the b^(th) combined heat and power unit in the integrated heat and electricity system at the t^(th) scheduling period, p_(x,t) represents an active power of an x^(th) thermal power unit in the integrated heat and electricity system at the t^(th) scheduling period, N_(TU) represents a number of thermal power units in a power grid in the integrated heat and electricity system, F_(TU)(p_(x,t)) represents a running cost of the x^(th) thermal power unit in the power grid in the integrated heat and electricity system, Δt represents an interval between two neighboring scheduling periods, and a value of Δt is 15 minutes.

(2) Equality constraints for a steady-state safe operation of the power grid and the heating network in the integrated heat and electricity system may be set, which may include the following.

(2-1) A power flow constraint of the power grid in the integrated heat and electricity system is denoted by formulas of:

${P_{i,t} = {U_{i,t}{\sum\limits_{j \in i}{U_{j,t}\left( {{G_{ij}{\cos \left( {\theta_{i,t} - \theta_{j,t}} \right)}} + {B_{ij}{\sin \left( {\theta_{t} - \theta_{j,t}} \right)}}} \right)}}}},i,{j = 1},2,{\cdots \mspace{14mu} n},{t = 1},2,\ldots \mspace{11mu},16$ ${Q_{i,t} = {U_{i,t}{\sum\limits_{j \in i}{U_{j,t}\left( {{G_{ij}{\sin \left( {\theta_{i,t} - \theta_{j,t}} \right)}} - {B_{ij}{\cos \left( {\theta_{i,t} - \theta_{j,t}} \right)}}} \right)}}}},i,{j = 1},2,{\cdots \mspace{14mu} n},{t = 1},2,\ldots \mspace{11mu},16,$

where, P_(i,t) represents an injected active power of a i^(th) node of the power grid in the integrated heat and electricity system at the t^(th) scheduling period,

_(i,t) represents an injected reactive power of the i^(th) node of the power grid in the integrated heat and electricity system at the t^(th) scheduling period, θ_(i,t) represents a voltage phase angle of the i^(th) node at the t^(th) scheduling period, θ_(j,t) represents a voltage phase angle of the j^(th) node at the t^(th) scheduling period, U_(i,t) represents a voltage amplitude of the i^(th) node at the t^(th) scheduling period, U_(j,t) represents a voltage amplitude of the j^(th) node at the t^(th) scheduling period, G_(ij) represents a real part of an element at i^(th) cow and j^(th) column in a node admittance matrix Y of the power grid, B_(ij) represents an imaginary part of the element at i^(th) cow and j^(th) column in the node admittance matrix Y of the power grid, the node admittance matrix Y of the power grid is obtained from an energy management system of the integrated heat and electricity system, and n represents a number of nodes of the power grid.

(2-2) A pipe pressure loss constraint of the heating network in the integrated heat and electricity system is denoted by a formula of:

ΔH _(l,t) =S _(l) m _(l,t) |m _(l,t) |, t=1,2, . . . ,16,

where, ΔH_(l,t) represents a pressure loss of a l^(th) pipe of the heating network in the integrated heat and electricity system at the t^(th) scheduling period, S_(l) represents a resistance characteristic coefficient, a range of S_(l) is 10 Pa/(kg/s)²≤S_(l)≤500 Pa/(kg/s)², and m_(l,t) represents a mass flow rate of the l^(th) pipe at the t^(th) scheduling period.

(2-3) A hydraulic characteristic constraint of a circulating pump of the heating network in the integrated heat and electricity system is denoted by a formula of:

H _(P,t) =H ₀ −S _(p) m _(p,t) ² t=1,2, . . . ,16,

where, H_(P,t) represents a head of delivery of the circulating pump at the t^(th) scheduling period, H₀ represents a static head of delivery of the circulating pump, S_(p) represents a resistance coefficient of the circulating pump, H₀ and S_(p) are obtained from a factory manual of the circulating pump, and m_(P,t) represents a mass flow rate that flowing through the circulating pump at the t^(th) scheduling period.

(2-4) A pipe heating loss constraint of the heating network in the integrated heat and electricity system is denoted by a formula of:

${T_{e,l,t} = {{\left( {T_{h,l,t} - T_{a,l,t}} \right)e^{- \frac{\lambda \; L_{l}}{{C_{p}m_{l}},t}}} + T_{a,l,t}}},{t = 1},2,\ldots \mspace{11mu},16,$

where, T_(e,l,t) represents a tail-end temperature of the l^(th) pipe of the heating network at the t^(th) scheduling period, T_(h,l,t) represents a head-end temperature of the l^(th) pipe at the t^(th) scheduling period, T_(a,l,t) represents an ambient temperature of the l^(th) pipe at the t^(th) scheduling period, m_(l,t) represents the mass flow rate of the l^(th) pipe at the t^(th) scheduling period, L₁ represents a length of the l^(th) pipe, C_(p) represents a specific heat capacity of water, a value of the specific heat capacity is 4182 Joules/(kg·° C.), λ represents a heat transfer coefficient per unit length of the pipe, λ is obtained from the energy management system of the integrated heat and electricity system, and e represents a natural logarithm.

(2-5) A temperature constraint of a junction node of a plurality of pipes of the heating network in the integrated heat and electricity system is denoted by a formula of:

${{\left( {\sum\limits_{l \in S_{n}^{out}}m_{l,t}} \right)T_{h,l,t}} = {{\sum\limits_{l \in S_{n}^{i\; n}}{m_{l,t}T_{e,l,t}}} - Q_{J,n,t}}},{t = 1},2,\ldots \mspace{11mu},16,$

where, m_(l,t) represents the mass flow rate of the l^(th) pipe at the t^(th) scheduling period, T_(e,l,t) represents the tail-end temperature of the l^(th) pipe of the heating network at the t^(th) scheduling period, T_(h,l,t) represents the head-end temperature of the l^(th) pipe at the t^(th) scheduling period, S_(n) ^(out) represents a set of all pipe numbers flowing out of a n^(th) node of the heating network, S_(n) ^(in) represents a set of all pipe numbers flowing into the n^(th) node of the heating network, and

_(j,n,t) represents a heating power of the n^(th) node of the heating network at the t^(th) scheduling period.

(2-6) A coupling constraint between the power grid and the heating network coupled by combined heat and power units in the integrated heat and electricity system is denoted by formulas of:

p _(b,t)=Σ_(k−1) ^(NK) ^(b) α_(b,t) ^(k) P _(b) ^(k) , q _(b,t)=Σ_(k=1) ^(NK) ^(b) α_(b,t) ^(k)

_(b) ^(k) , t=1,2, . . . ,16,

where, p_(b,t) represents the active power of the b^(th) combined heat and power unit at the t^(th) scheduling period, q_(b,t) represents the heating power of the b^(th) combined heat and power unit at the t^(th) scheduling period, (P_(b) ^(k),

_(b) ^(k)) represents the k^(th) vertex of an approximate convex polygon in a running feasible region of the b^(th) combined heat and power unit, α_(b,t) ^(k) represents a k^(th) combination coefficient of the b^(th) combined heat and power unit at the t^(th) scheduling period, Σ_(k−1) ^(NK) ^(b) α_(b,t) ^(k)=1, 0≤α_(b,t) ^(k)1, NK_(b) represents a number of vertexes of the approximate polygon in the running feasible region of the b^(th) combined heat and power unit, and the approximate polygon in the running feasible region of the combined heat and power unit is obtained from a factory instruction of the combined heat and power unit.

(2-7) A coupling constraint between the power grid and the heating network coupled by a circulating pump in the integrated heat and electricity system is denoted by a formula of:

${P_{P,t} = \frac{m_{P,t}gH_{P,t}}{10^{6}\eta_{P}}},{t = 1},2,\ldots \mspace{11mu},16,$

where, P_(P,t) represents an active power consumed by the circulating pump at the t^(th) scheduling period, g represents an acceleration of gravity, η_(P) represents an efficiency of the circulating pump, a value of η_(P) ranges from 0 to 1, m_(P,t) represents the mass flow rate that flowing through the circulating pump at the t^(th) scheduling period, and H_(P,t) represents the head of delivery of the circulating pump at the t^(th) scheduling period.

(2-8) A coupling constraint between the power grid and the heating network coupled by a heat pump in the integrated heat and electricity system is denoted by a formula of:

P _(hp,t) =C _(hp)

_(hp,t) ,t=1,2, . . . ,16,

where,

_(hp,t) represents a heating power emitted by the heat pump in the integrated heat and electricity system at the t^(th) scheduling period, p_(hp,t) represents an electric power consumed by the heat pump at the t^(th) scheduling period, C_(hp) represents a heat generation efficiency of the heat pump, and C_(hp) is obtained from a factory manual of the heat pump.

(3) Inequality constraints for a steady-state safe operation of the power grid and a heating network in the integrated heat and electricity system may be set, which may include the following.

(3-1) The voltage amplitude U_(i,t) of the i^(th) node of the power grid in the integrated heat and electricity system at the t^(th) scheduling period is between a set upper and lower limits U _(i), Ū_(i) of safe operating voltage of the power grid, in which, U _(i) is 0.95 times a rated voltage of the i^(th) node, and Ū_(i) is 1.05 times the rated voltage of the i^(th) node, i.e.:

U _(i) ≤U _(i,t) ≤Ū _(i) ,t=1,2, . . . ,16.

(3-2) A transmission capacity S_(l,t) of a l^(th) branch of the power grid in the integrated heat and electricity system at the t^(th) scheduling period is less than or equal to a set maximum value S _(l) of the transmission capacity for safe operating of the power grid, i.e.,:

S _(l,t) ≤S _(l) , t=1,2, . . . ,16.

(3-3) Ramp constraints of active powers of combined heat and power units of the power grid of the integrated heat and electricity system are denoted by a formula of:

−Δp _(b) ^(down) ·Δt≤p _(b,t) −p _(b,t−1) ≤Δp _(b) ^(up) ·Δt,t=1,2, . . . ,16,

where, Δp_(b) ^(up) represents an upward ramp rate of the active power of the b^(th) combined heat and power unit, Δp_(b) ^(down) represents a downward ramp rate of the active power of the b^(th) combined heat and power unit, Δp_(b) ^(up) and Δp_(b) ^(down) are obtained from a factory instruction of the combined heat and power unit, Δt represents the interval between two neighboring scheduling periods, p_(b,t) represents the active power of the b^(th) combined heat and power unit at the t^(th) scheduling period, and p_(b,t−1) represents the active power of the b^(th) combined heat and power unit at the (t−1)^(th) scheduling period.

(3-4) Ramp constraints of heating powers of combined heat and power units of the power grid of the integrated heat and electricity system are denoted by a formula of:

−Δq_(d) ^(down) ·Δt≤q _(b,t) −q _(b,t−1) ≤Δq _(b) ^(up) ·Δt,t=1, 2, . . . ,16,

where, Δq_(b) ^(up) represents an upward ramp rate of the heating power of the b^(th) combined heat and power unit,Δq_(b) ^(down) represents a downward ramp rate of the heating power of the b^(th) combined heat and power unit, Δq_(b) ^(up) and Δq_(b) ^(down) are obtained from a factory instruction of the combined heat and power unit, Δt represents the interval between two neighboring scheduling periods, q_(b,t) represents the heating power of the b^(th) combined heat and power unit at the t^(th) scheduling period, and q_(b,t−1) represents the heating power of the b^(th) combined heat and power unit at the (t−1)^(th) scheduling period.

(3-5) Ramp constraints of active powers of thermal power units of the power grid of the integrated heat and electricity system are denoted by a formula of:

Δ−p _(x) ^(down) ·Δt≤p _(x,t) −p _(x,t−1) ≤Δp _(x) ^(up) ·Δt,t=1,2, . . . ,16,

where, Δp_(x) ^(up) represents an upward ramp rate of the active power of the x^(th) thermal power unit, Δp_(x) ^(down) represents a downward ramp rate of the active power of the x^(th) thermal power unit, Δp_(x) ^(up) and Δp,_(x) ^(down) are obtained from a factory instruction of the thermal power unit, Δt represents the interval between two neighboring scheduling periods, p_(x,t) represents the active power of the x^(th) thermal power unit at the t^(th) scheduling period, and p_(x,t−1) represents the active power of the x^(th) thermal power unit at the (t−1)^(th) scheduling period.

(3-6) The active power p_(b,t) of the b^(th) combined heat and power unit at the t^(th) scheduling period is between a set upper and lower limits P _(b) and P _(b) of the active power of the b^(th) combined heat and power unit, i.e.,

P _(b) ≤P _(b,t) ≤p _(b) ,t=1,2, . . . ,16 .

(3-7) The heating power q_(b,t) of the b^(th) combined heat and power unit at the t^(th) scheduling period is between a set upper and lower limits q _(b) and q _(b) of the heating power of the b^(th) combined heat and power unit, i.e.,

p _(b) ≤p _(b,t) ≤p _(b) ,t=1,2, . . . ,16.

(3-8) The active power p_(x,t) of the x^(th) thermal power unit at the t^(th) scheduling period is between a set upper and lower limits p _(x), and p _(x) of the active power of the x^(th) thermal power unit, i.e.,

p _(x) ≤p _(x,t) ≤p _(x) ,t=1,2, . . . ,16.

(3-9) The mass flow rate m_(l,t) of the l^(th) pipe at the t^(th) scheduling period is less than or equal to an upper limit m _(l) of the mass flow rate for safe operation of the heating network, i.e.:

T≤T _(i,t) ≤T,t=1,2, . . . ,16.

(3-10) A return water temperature T_(i,t) of an i^(th) heat exchange station in the heating network at the t^(th) scheduling period is between a set upper and lower limits T and T of the return water temperature for safe operation of the heating network, i.e.:

T≤T _(i,j) T,t=1,2, . . . , 16.

(4) Using an interior point method, the equation in the above step (1) is used as the objective function, and all equations in the above step (2) and step (3) are used as the constraints, to solve to obtain the active power and the heating power of each combined heat and power unit, the active power of each thermal power unit, the heating power of each heat pump, and the active power consumed by each circulating pump, as the intra-day rolling scheduling scheme of the integrated heat and electricity system, realizing the intra-day rolling scheduling of the integrated heat and electricity system.

The intra-day rolling scheduling method for the integrated heat and electricity system, proposed by the present disclosure, has the characteristics and effects as follows. The method considers the mutual influence of the integrated heat and electricity system and realizes the intra-day rolling scheduling method for the integrated heat and electricity system. This method can continuously modify the current intra-day rolling scheduling plan to obtain a more accurate intra-day rolling scheduling scheme. Compared with independent optimization analysis of the power grid and the heating system, it may get a better scheduling scheme (lower total operating cost), and improve scheduling flexibility. The method can be applied to the formulation of the intra-day rolling scheduling method for the integrated heat and electricity system, which is beneficial to improving the energy efficiency of the integrated heat and electricity system and reducing operating costs. 

What is claimed is:
 1. An intra-day rolling scheduling method for an integrated heat and electricity system, comprising: establishing an objective function for scheduling of the integrated heat and electricity system, the objective function aiming to make operating costs of the integrated heat and electricity system to be a minimum, the integrated heat and electricity system comprising combined heat and power units, thermal power units, heat pumps, and circulating pumps; establishing constraints for a steady-state safe operation of the integrated heat and electricity system; and solving the objective function based on the constraints by an interior point method, to obtain an active power and a heating power of each combined heat and power unit, an active power of each thermal power unit, a heating power of each heat pump, and an active power consumed by each circulating pump, as an intra-day rolling scheduling scheme of the integrated heat and electricity system.
 2. The method of claim 1, wherein the objective function for scheduling of the integrated heat and electricity system is denoted by a formula of: ${\min {\sum\limits_{t = 1}^{16}{\left( {{\sum\limits_{b = 1}^{N}{F\ \left( {p_{b,t},\ q_{b,t}} \right)}} + {\sum\limits_{x = 1}^{N_{TU}}{F_{TU}\left( p_{x,t} \right)}}} \right)\Delta t}}},$ where, p_(b,t) represents an active power of a b^(th) combined heat and power unit in the integrated heat and electricity system at a t^(th) scheduling period, q_(b,j) represents a heating power of the b^(th) combined heat and power unit at the t^(th) scheduling period, N represents a number of combined heat and power units, F(p_(b,t),q_(b,t)) represents a running cost of the b^(th) combined heat and power unit at the t^(th) scheduling period, p_(x,t) represents an active power of an x^(th) thermal power unit in a power grid in the integrated heat and electricity system at the t^(th) scheduling period, N_(TU) represents a number of thermal power units, F_(TU)(p_(x,t)) represents a running cost of the x^(th) thermal power unit, and Δt represents an interval between two neighboring scheduling periods.
 3. The method of claim 2, wherein the constraints comprise: a power flow constraint of the power grid, a pipe pressure loss constraint of a heating network in the integrated heat and electricity system, a hydraulic characteristic constraint of a circulating pump of the heating network, a pipe heating loss constraint of the heating network, a temperature constraint of a junction node of a plurality of pipes of the heating network, a coupling constraint between the power grid and the heating network coupled by combined heat and power units, a coupling constraint between the power grid and the heating network coupled by the circulating pump, a coupling constraint between the power grid and the heating network coupled by the heat pump, and upper and lower limit constraints.
 4. The method of claim 3, wherein the power flow constraint of the power grid is denoted by formulas of: ${P_{i,t} = {U_{i,t}{\sum\limits_{j \in i}{U_{j,t}\left( {{G_{ij}{\cos \left( {\theta_{i,t} - \theta_{j,t}} \right)}} + {B_{ij}{\sin \left( {\theta_{t} - \theta_{j,t}} \right)}}} \right)}}}},i,{j = 1},2,{\cdots \mspace{14mu} n},{t = 1},2,\ldots \mspace{11mu},16$ ${Q_{i,t} = {U_{i,t}{\sum\limits_{j \in i}{U_{j,t}\left( {{G_{ij}{\sin \left( {\theta_{i,t} - \theta_{j,t}} \right)}} - {B_{ij}{\cos \left( {\theta_{i,t} - \theta_{j,t}} \right)}}} \right)}}}},i,{j = 1},2,{\cdots \mspace{14mu} n},{t = 1},2,\ldots \mspace{11mu},16,$ where, P_(i,t) represents an injected active power of a i^(th) node of the power grid at the t^(th) scheduling period,

_(l,t) represents an injected reactive power of the i^(th) node of the power grid at the t^(th) scheduling period, θ_(i,t) represents a voltage phase angle of the i^(th) node at the t^(th) scheduling period, θ_(j,t) represents a voltage phase angle of the j^(th) node at the t^(th) scheduling period, U_(i,t) represents a voltage amplitude of the i^(th) node at the t^(th) scheduling period, U_(j,t) represents a voltage amplitude of the j^(th) node at the t^(th) scheduling period, G represents a real part of an element at i^(th) cow and j^(th) column in a node admittance matrix Y of the power grid, B_(ij) represents an imaginary part of the element at i^(th) cow and j^(th) column in the node admittance matrix Y of the power grid, and n represents a number of nodes of the power grid.
 5. The method of claim 3, wherein the pipe pressure loss constraint of the heating network is denoted by a formula of: ΔH _(l,t) =S _(l) m _(l,t) |m _(l,t) |,t=1,2, . . . ,16, where, ΔH_(i,t) represents a pressure loss of a l^(th) pipe of the heating network at the t^(th) scheduling period, S_(l) represents a resistance characteristic coefficient, and m_(i,t) represents a mass flow rate of the l^(th) pipe at the t^(th) scheduling period.
 6. The method of claim 3, wherein the hydraulic characteristic constraint of a circulating pump of the heating network is denoted by a formula of: H _(P,t) =H ₀ −S _(p) m _(p,t) ² ,t=1,2, . . . ,16, where, H_(p,t) represents a head of delivery of the circulating pump at the t^(th) scheduling period, H₀ represents a static head of delivery of the circulating pump, S_(p) represents a resistance coefficient of the circulating pump, and m_(P,t) represents a mass flow rate that flowing through the circulating pump at the t^(th) scheduling period.
 7. The method of claim 3, wherein the pipe heating loss constraint of the heating network is denoted by a formula of: ${T_{e,l,t} = {{\left( {T_{h,l,t} - T_{a,l,t}} \right)e^{- \frac{\lambda \; L_{l}}{{C_{p}m_{l}},t}}} + T_{a,l,t}}},{t = 1},2,\ldots \mspace{11mu},16,$ where, T_(e,l,t) represents a tail-end temperature of the l^(th) pipe of the heating network at the t^(th) scheduling period, T_(h,l,t) represents a head-end temperature of the l^(th) pipe at the t^(th) scheduling period, T_(a,l,t) represents an ambient temperature of the l^(th) pipe at the t^(th) scheduling period, m_(l,t) represents the mass flow rate of the l^(th) pipe at the t^(th) scheduling period, L_(l) represents a length of the l^(th) pipe, C_(p) represents a specific heat capacity of water, λ represents a heat transfer coefficient per unit length of the pipe, and e represents a natural logarithm.
 8. The method of claim 3, wherein the temperature constraint of the junction node of the plurality of pipes of the heating network is denoted by a formula of: ${{\left( {\sum\limits_{l \in S_{n}^{out}}m_{l,t}} \right)T_{h,l,t}} = {{\sum\limits_{l \in S_{n}^{i\; n}}{m_{l,t}T_{e,l,t}}} - Q_{J,n,t}}},{t = 1},2,\ldots \mspace{11mu},16,$ where, m_(l,t) represents a mass flow rate of the l^(th) pipe at the t^(th) scheduling period, T_(e,l,t) represents a tail-end temperature of the l^(th) pipe of the heating network at the t^(th) scheduling period, T_(h,l,t) represents a head-end temperature of the l^(th) pipe at the t^(th) scheduling period, S_(n) ^(out) represents a set of all pipe numbers flowing out of a n^(th) node of the heating network, S_(n) ^(in) represents a set of all pipe numbers flowing into the n^(th) node of the heating network, and

_(j,n,t) represents a heating power of the n^(th) node of the heating network at the t^(th) scheduling period.
 9. The method of claim 3, wherein the coupling constraint between the power grid and the heating network coupled by the combined heat and power units is denoted by formulas of: p _(b,t)Σ_(k=1) ^(NK) ^(b) α_(b,t) ^(k) P _(b) ^(k) , q _(b,t)=Σ_(k=1) ^(NK) ^(b) α_(b,t) ^(k)

_(b) ^(k) , t=1,2, . . . ,16, where, (P_(b) ^(k),

_(b) ^(k)) represents a k^(th) vertex of an approximate convex polygon in a running feasible region of the b^(th) combined heat and power unit, α_(b,t) ^(k) represents a k^(th) combination coefficient of the b^(th) combined heat and power unit at the t^(th) scheduling period, Σ_(k=1) ^(NK) ^(b) α_(b,t) ^(k)=1, 0≤α_(b,t) ^(k)≤1, and NK_(b) represents a number of vertexes of the approximate polygon in the running feasible region of the b^(th) combined heat and power unit.
 10. The method of claim 3, wherein the coupling constraint between the power grid and the heating network coupled by the circulating pump is denoted by a formula of: ${P_{P,t} = \frac{m_{P,t}gH_{P,t}}{10^{6}\eta_{P}}},{t = 1},2,\ldots \mspace{11mu},16,$ where, P_(P,t) represents an active power consumed by the circulating pump at the t^(th) scheduling period, g represents an acceleration of gravity, η_(P) represents an efficiency of the circulating pump, represents a mass flow rate that flowing through the circulating pump at the t^(th) scheduling period, and H_(P,t) represents a head of delivery of the circulating pump at the t^(th) scheduling period.
 11. The method of claim 3, wherein the coupling constraint between the power grid and the heating network coupled by the heat pump is denoted by a formula of: P _(hp,t) =C _(hp)

_(hp,t) ,t=1,2, . . . ,16, where, Q_(hp,t) represents a heating power emitted by the heat pump at the t^(th) scheduling period, p_(hp,t) represents an electric power consumed by the heat pump at the t^(th) scheduling period, and C_(hp) represents a heat generation efficiency of the heat pump.
 12. The method of claim 4, wherein upper and lower limit constraints comprise: U _(i)≤U_(i,t)≤Ū_(i),t=1,2, . . . ,16, where, U represents an upper limit for safe operating voltage of the i^(th) node, and Ū_(l) represents a lower limit for safe operating voltage of the i^(th) node; S_(l,t)≤S _(l),t=1,2, . . . ,16, where, S_(l,t) represents a transmission capacity of a l^(th) branch of the power grid, and S _(l) represents a maximum transmission capacity of the power grid; −Δp_(b) ^(down)·Δt≤p_(b,t)−p_(b,t−1)≤Δp_(b) ^(up)·Δt,t=1, 2, . . . ,16 , where, Δp_(b) ^(up) represents an upward ramp rate of the active power of the b^(th) combined heat and power unit, Δp_(b) ^(down) represents a downward ramp rate of the active power of the b^(th) combined heat and power unit, and P_(b,t−1) represents an active power of the b^(th) combined heat and power unit at the (t−1)^(th) scheduling period; −Δq_(b) ^(down)·Δt≤≤q_(b,t)−q_(b,t−1)≤Δq_(b) ^(up)·Δt,t=1,2, . . . ,16, where, Δq_(b) ^(up) represents an upward ramp rate of the heating power of the b^(th) combined heat and power unit, Δq_(b) ^(up) represents a downward ramp rate of the heating power of the b^(th) combined heat and power unit, and q_(b,t−1) represents a heating power of the b^(th) combined heat and power unit at the (t−1)^(th) scheduling period; −Δp_(x) ^(down)·Δt≤p_(x,t)−p_(x,t=1)≤Δp_(x) ^(up)·Δt,t=1,2, . . . ,16, where, Δp_(x) ^(up) represents an upward ramp rate of the active power of the x^(th) thermal power unit, Δp_(x) ^(down) represents a downward ramp rate of the active power of the x^(th) thermal power unit, and p_(x,t−1) represents an active power of the x^(th) thermal power unit at the (t−1)^(th) scheduling period; p _(b)≤p_(b,t) P _(b),t=1,2, . . . ,16, where, p _(b) represents an upper limit of the active power of the b^(th) combined heat and power unit and p _(b) represents a lower limit of the active power of the b^(th) combined heat and power unit; q _(b) ≤q _(b,t) ≤q _(b) ,t=1,2, . . . ,16, where, q _(b) represents an upper limit of the heating power of the b^(th) combined heat and power unit and q _(b) represents a lower limit of the heating power of the b^(th) combined heat and power unit; and p _(x)≤p_(x,t)≤p _(x),t=1,2, . . . ,16, wherein, p _(x) represents an upper limit of the active power of the x^(th) thermal power unit and p _(x) represents a lower limit of the active power of the x^(th) thermal power unit.
 13. The method of claim 5, wherein upper and lower limit constraints comprise: 0≤m_(l,t)≤m _(l),t=1,2, . . . ,16, where, m _(l) represents an upper limit of the mass flow rate of the heating network; and T≤T_(i,j≤T,t=)1,2, . . . ,16, where, T_(i,j) represents a return water temperature of an i^(th) heat exchange station in the heating network at the t^(th) scheduling period, T represents an upper limit of the return water temperature of the heating network, and T represents a lower limit of the return water temperature of the heating network. 